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Abstract 

Fast Fourier Transform (FFT) is an efficient algorithm to compute the Discrete Fourier 
Transform (DFT) and its inverse. In this paper, we pay special attention to the description 
of complex-data FFT. We analyze two common descriptions of FFT and propose a new 
presentation. Our heuristic description is helpful for students and programmers to grasp 
the algorithm entirely and deeply. 
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1 Introduction 

Denote the set of all complex numbers by C. Let oj = e 2/Kt l n be a primitive nth root of unity. 
Let 

/(*)= Y. fij£C[x] (1) 

0<i<n 

be a polynomial of degree less than n with its coefficient vector (/o, • • • , f n -i) £ C n - The 
C-linear map 



DFT W : 



C n _^ C n 



which evaluates a polynomial at the powers of co is called the Discrete Fourier Transform. 
Apparently, the Discrete Fourier Transform is a special multipoint evaluation at the powers 

1,UJ,L0 2 , • • • ,W n_1 . 

A beginner will often evaluate the polynomial (1) in a manner that: first / n _ix n_1 is 
calculated, then f n -2X n ~ 2 , • • • , fix, and finally all of the terms of (1) are added together. 
Such a process involves lots of multiplications and additions. There is an elegant way to 
evaluate a polynomial, which is called Horner's rule. It evaluates f(x) as follows: Rearrange 



this computation as 

f{x) = {■■■ (f n - lX + / n _ 2 )x + ..-)x + f (2) 

Then start with / n _i, multiply by x, add f n -2, multiply by x, multiply by x, add /o- 
If we apply the Horner's rule to compute DFT W , it needs to run the evaluation algorithm n 
times. Can we find more efficient algorithm to compute Discrete Fourier transform? 

Fast Fourier Transform (FFT) is an efficient algorithm to compute the Discrete Fourier 
Transform (DFT) and its inverse. It was popularized by a publication of J. Cooley and 
J. Tukey PQ in 1965. It has been called the most important numerical algorithm of our 
lifetime. In this paper, we will pay special attention to the description of complex-data FFT. 
We will analyze two common descriptions of FFT [3 0] and propose an explicit and heuristic 
presentation of FFT. 



2 Analysis of two common descriptions of FFT 
2.1 Description-1 

We refer to Ref.[3] for the first description of FFT. 

Let n = 2 k £ M with k 6 j\f, to 6 C be a primitive nth root of unity, and f(x) = 
J2o<i<nfi xl ^ of degree less than n. To evaluate f(x) at the powers l,w,w 2 , • • • ,w n_1 , 
we divide fix) by x n l 2 — 1 and x n l 2 + 1 with remainder: 

f(x) = q {x){x n l 2 - 1) + r (x) = qi {x){x n ' 2 + 1) + n (x) (3) 

for some qo(x), ro(x), qi(x), r\{x) 6 C[x] of degree less than n/2. 
Plugging in a power of u for x in (3), we find 

f(co 2e+1 ) = qi {uo 2l+1 ){uj nt uj n/2 + 1) + ri(w M+1 ) = ri (u 2e+1 ) 

for all < £ < n/2. We have used the facts that uj n = 1 and uf 1 ! 2 = — 1. It remains to evaluate 
vq(x) at the even powers of u and r±(x) at the odd powers. Now uj 2 is a primitive (n/2)th 
root of unity, and hence the first task is a DFT of order n/2. But also the evaluation of r\[x) 
can be reduced to a DFT of order n/2 by noting that ri(uj 2 ^ +l ) = r\{oj 2t ) for r\{x) = r\{ujx). 
Since n is a power of 2, we can proceed recursively to evaluate vq(x) and r\{x) at the powers 
1, uj 2 , ■ ■ ■ , ui 2n ~ 2 of to 2 , and obtain the following algorithm. 
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Table 1: Description- 1 



Input: n = 2 k £ N with k £ M, f(x) = J2o<i<n fi x ' 1 ^ C[x], and 

the powers lj, uj 2 , • • • , of a primitive nth root of unity to G C. 
Output: DFT W (/) = (/(l), /(w), • • • , /(w"" 1 )) G C. 

1. if n = 1 then return /o 

2. r (s) « E <j<n/2(/j + fj+n/2)x J , rj(x) i E <j<n/2(/i ~ /j+n/2 V X 3 

3. call the algorithm recursively to evaluate vq(x) and rj(x) at the powers of u 2 

4. return (r (l), r (u; 2 ), rj(w 2 ), • • • , r (cj n ~ 2 ), r\(uj n ~ 2 )) 



2.2 Analysis of Description-1 

To investigate the working flow in description-1, we take the following polynomial as an 
example: 

f(x) = f + fix + f 2 x 2 + --- + f 7 x 7 €C[x], co = e 2 ^ 8 (4) 
At first, we have the following two remainders: 

ro(x) = (k + k) + (fi + k)x + (k + k)x 2 + (k + k)x 3 (5) 
rl(x) = (f Q -f 4 ) + (f 1 -f 5 ) u >x + (f 2 -f 6 )u 2 x 2 + (f 3 -f 7 )cu 3 x 3 (6) 

If we evaluate them at (uj 2 ) , then r (l) = /(l), r*(l) = f(u). 

We now proceed to the step 3. By (5) we have the following two remainders: 

[(/o + k) + (k + /e)] + [(/i + /s) + (/s + / 7 )]x (7) 
[(/o + U) - (h + k)] + [(/i + /s) - (/s + fi)]u x (8) 

Likewise, by (6) we have the following two remainders: 

[(/o - U) + (/ 2 - k)u 2 ] + [(/i - / 5 )w + (/ 3 - / 7 )w 3 ]x (9) 
[(/o - /4) - (k - k)co 2 ] + [(/i - / 5 )w - (/ 3 - M^ 3 ]^ (10) 

Should we evaluate them at (w 2 ) 1 ? If that we can not find correct answers. 

Drawback: The sentence that, "call the algorithm recursively to evaluate ro(x) and r\{x) 
at the powers of w 2 " , is really too vague to specify the working flow. 
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2.3 Description-2 



We refer to Ref.[4] for the second description of FFT. 

Let n = 2 k £ N with k £ M, to € C be a primitive nth root of unity, and f(x) = 
So<i<n/« x * ^ ^t x ] °f degree less than n. To evaluate /(x) at the powers 1,oj,w 2 , • • • ,u n ~ 1 , 
we carry out the following scheme. (In these formulas the parameters Sj and tj are either 
or 1, so that each "pass" represents 2 k elementary computations.) 

Table 2: Description-2 

Pass 0. Let A^\t k _i, ■ ■ ■ ,t ) = ft, where t = (t k -i, ■ ■ ■ ,t ) 2 . 
Pass 1. Set AW(s k -i,t k -2,---,to) <- 

A® (0, t fc _ 2) • • • , t ) + v 2 "- 1 **-! A® (l,t k - 2 , ■ ■ ■ , t ) 
Pass 2. Set A® (s k -i, s k _ 2 , ifc-3, • • • , t ) +- 

AW(«fc-l, 0, fjfc-3, • • • , to) + c^-^-^-^AH^, 1, t fc _ 3 , • • • , t ) 

Pass k. Set A^ (s k -i, ■ ■ ■ , si , s ) ^~ 

A^ k ~^{s k -i, ■ ■ • , si, 0) + w («o'i"-k-i^[*-i]( gw , ■ ■ ■ , ai, 1) 

It is easy to find that 

A^(s k -i, ■ ■ ■ ,si,s ) = f{u s ), where s = (s , si, • • • , s fc _i) 2 (11) 
Notice that the binary digits of s are reversed in the final result (11). 



2.4 Analysis of Description-2 

To investigate the working flow in description-2, we also take the polynomial (4) as an example. 
Clearly, we have 

A[°](000) = / o , Al°](001) = / 1 , Al°](010) = / 2 , ^(011) =/ 3 
A[°](100) = / 4 , ^°](101)=/ B) A [0] (H0) - J 6 , ^°](lll)=/ 7 



In pass 1, 



AM(OOO) = A™ (000) +^[°1(100) = fo + U, AW(100) = A™ (000) + oj 2 ' A® (100) = f - f 4 

AM(OOl) = A™ (001) + A^lOl) = A + / B , AW(101) = A™ (001) + w 2 ' A™ (101) = A - / 5 

AW(010) = A[°l(010) + A[°l(110) = h + fa, AW(llO) = A™ (010) + cu 2 ' A® (110) = f 2 - f e 

A™ (011) = AM (011)+ AM (111) = h + f 7 , AW(lll) = A[°](011)+ w 22 A[°](ill) =/ 3 -/ 7 
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In pass 2, 



In pass 3, 



,4 1 2 1 (000) = AW(000) + w 2x(00)2 AW(010) 
^(OOl) = AW(OOl) +w 2x (° 0)2 AW(011) 
A^OIO) = AW(OOO) + w 2x ( 10 ' 2 A[ : 
^(Oll) = AW(001)+w 2x(10)2 
A^lOO) = A[ 1 ](100)+w 2x ( 01 ' 2 

^(ioi) = aW(ioi) + w 2x ( 01 ) 2 

A^llO) = AW(100)+w 2x ( n ) 2 

aP](ih) ^^(ioij + w 2 ** 11 ^ 



(010) 
A [1] (011) 
AW(llO) 
AW(lll) 
AW(llO) 

(111) 



(/o + h) + (h + fe) 

(.h + h) + (h + h) 

(/O + k) - C/2 + h) 

(fi + h) - (h + h) 
(/o- U) + u\f2- h) 

(h- h)+u 2 {h- h) 

(/o-/4)-w 2 (/ 2 -/ 6 ) 
(fi-h)-^(fs-fr) 



^(OOO) 


= A 121 (000) + 


W (000) M [2] (OQ1) = 


[(/o 


+ / 4 ) + 


(/ 2 + /e)] + [(A + h 


) + {h + h)] = 


/(i) 


^(OOl) 


= A [21 (000) + 


W (100) M [2] (001) = 


K/o 


+ U) + 


(/2 + /e)]-[(/i+/ 5 


) + {h + h)] = 


/(^ 4 ) 


a[ 3 ](oio) 


= A 121 (010) + 


w (010) M [2] (on) = 


K/o 


+ U)- 


(/ 2 + / 6 )]+w 2 [(/i + 


h) - (h + h)] 


= /(- 2 ) 


^(Oll) 


= ,4 [21 (010) + 


w (iio) 2A [2] (on) = 


K/o 


+ h)- 


(/ 2 + / 6 )]-w 2 [(/i + 


h) - (h + h)\ 


= /(^ 6 ) 


^[31(100) 


= A 121 (100) + 


w (ooi) M [2] (101) = 


K/o 


-h) + 


^ 2 (/ 2 -/ 6 )]+^[(/i- 


-h) + u 2 (h- 


/t)] = /M 


^(lOl) 


= A [21 (100) + 


w (ioi) 2j4 [2] (101) = 


[(/o 


-h) + 


^(/ 2 -/6)]-^[(/l- 


-h) + u 2 (h- 


/r)] = Z(" 5 ) 


^(llO) 


= ,4[ 2 1(110) + 


w (oii) M [2] (ln) = 


[(/o 


-k)- 


^ 2 (/ 2 -/6)]+^[(/l 




-/ 7 )] = /(^ 3 ) 


^(lll) 


= Al 2 l(110) + 


w (lll)2 A [2] (m) = 


[(/o 


-h)- 


^ 2 (/ 2 -/6)]-^ 3 [(/l 


- f5 )-u 2 (f s - 


-/ 7 )] = /(^) 



Drawback: The description-2 almost veils the elegant idea behind FFT. 



3 A heuristic description of FFT 

We now present an explicit and heuristic description of FFT, which can be regarded as the 
combination of the above two descriptions. 

The basic idea of the new description is to recursively split polynomials and shrink coef- 
ficients. Concretely, given a polynomial f(x) = / + fix + /2a; 2 + • • • + /7a; 7 £ C[x] and a 
primitive 8th root of unity u = e 2ni/8 , we split it into two polynomials of degree less than 4. 
Their coefficients can be individually obtained by shrinking the symmetric two coefficients of 
original polynomial. See the following working flow for details. 



5 




6 



Refining the method, we have the following heuristic description for FFT: 

Table 3: Description-3 



Input: n = 2 k G Af with k ejV, f(x) = Eo<i<n fi x% and 

the powers uj, uj 2 , • • • , w n_1 of a primitive nth root of unity wGC. 
Output: DFT W (/) = (/(l), /(w), • • • , /(w"" 1 )) G C*. 

Pass 1. Set = ( E °^ <n/2 ^' + W | ^ . &1 = ° 

1 E <i<n/2 (/J' - fj+n/2) (w^, &1 = 1 

Pass 2. For each 6j G {0, 1}, 1 < i < 2, compute 

Eo<K^(/] 6ll +/]+ 1 I/ 22 )^' ^ = 

Eo^n/^^-Zj+V)^ 2 ^' 62 = 1 



J[6i6 2 ]( 



ife 2 ---&fe] 



(x) 



Pass k. For each hi G {0, 1}, 1 < i < fc, compute 

/ f [6i62-6fc-i] 1 /•[&i&2---&fc-i]\ j 7 _ n 

Eo<i<n/2 fc (Jj + i i+n /2fe J X ' fe - U 

/ f [6i62-6fc-i] f [&i&2-6fc-i]\ / J 2 k - 1 „\j l _ 1 

L0<j<n/2* (/j ~ Jj+n/2>° J ^ X J' bk ~ i 

Reversal permutation. For each 6j G {0, 1}, 1 < i < k, set 



It is easy to prove the correctness of above method by induction and show that the 
algorithm requires kn complex-number additions. Incidentally, we refer to Ref.[2] for an 
explicit pseudocode of FFT (ITERATIVE-FFT) 

4 Conclusion 

In this paper, we analyze two common descriptions of FFT and propose a heuristic pre- 
sentation of complex-data FFT. We think the new presentation is helpful for students and 
programmers to grasp the method entirely and deeply. 
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